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1.  Introduction 


Metallic  fragments  occasionally  are  found  imbedded  in  combustible  cartridge  cases  as 
a  result  of  manufacturing  anomalies.  During  the  interior-ballistic  cycle,  these  particles  are 
released  into  the  interior  flow  field  as  the  case  is  consumed.  Components  of  the  propelling 
charge,  such  as  igniters,  may  also  be  a  source  of  metallic  fragments  during  normal 
functioning.  A  concern  of  system  designers  is  that  these  fragments  may  not  be  ejected  from 
the  gun  or  vaporized  and  pose  an  ignition  threat  either  to  subsequently  loaded  cases  or  to 
fuel-rich  combustion  gases  mixing  with  air  upon  opening  of  the  breech.  This  report  is  a  first 
effort  to  quantify  the  temperature  history  of  such  particles  in  a  120  mm  tank  cannon  as  a 
function  of  their  mass,  shape,  composition,  and  initial  location.  It  should  be  noted  that, 
while  the  fragment  temperature  history  is  prerequisite  to  establishing  the  seriousness  of  the 
ignition  threat  for  a  given  set  of  circumstances,  experiments  to  establish  the  ignitibility  of 
combustible  case  material  or  combustion  gases  are  required  to  complete  the  assessment. 

The  objective  of  the  work  documented  in  this  report  was  to  formulate  the  problem 
mathematically  and  to  solve  it  using  a  level  of  approximation  expedient  to  gaining  an 
appreciation  of  the  various  factors  involved,  given  a  one-month  time  constraint.  The  analysis 
shows  that  the  fragment  temperature  is  critically  dependent  upon  the  description  of  the 
particle  drag,  which  for  a  tumbling  wafer-like  geometry  is  difficult  to  assess.  This  uncertainty 
is  handled  by  computing  upper  and  lower  bounds  on  the  drag  and  determining  how  these 
bounds  map  into  the  solution  for  the  temperature  history.  Analytical  approximations  to  the 
equations  are  utilized,  and  no  doubt  the  temperature  histories  could  be  improved  by  solving 
the  equations  numerically.  It  is  believed,  however,  that  the  uncertainties  in  the  particle  drag 
overshadow  inaccuracies  introduced  by  the  mathematical  approximations. 


2.  Basic  Assumptions 

1.  Fragments  are  thin  wafers  such  that  the  surface  area  associated  with  their  edges  is 
negligible  compared  with  that  of  their  faces.  This  geometry  is  expected  and  simplifies  the 
mathematics  but  the  assumption  could  readily  be  removed  if  desired.  In  addition,  the 
fragments  are  assumed  to  be  square  wafers  for  convenience. 

2.  The  particles  are  not  in  thermal  contact  with  the  gun-tube  wall.  Presumably,  a  small 
fragment  would  be  kept  in  motion  by  turbulence  and  would  experience  only  brief  and 
occasional  contact  with  the  wall. 

3.  Time  for  intra-fragment  thermal  relaxation  is  short  compared  to  the  rate  of  heat  transfer 
to  the  fragment  (due  to  the  high  thermal  diffusivity  of  metal);  therefore,  no  temperature 
gradients  exist  in  the  fragment. 

4.  Radiation  heat  transfer  between  the  gas  and  the  particles  is  negligible.  (This  is  shown  to 
be  a  valid  approximation.) 
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5.  The  interior-ballistic  flow  is  adequately  characterized  by  the  IBHVG2  code.  Since  the 
time  interval  of  interest  is  after  flamespreading  is  complete,  this  lumped-parameter  code 
should  give  a  reasonable  first-order  description  of  the  flow.  An  assumption  which  is 
uncertain,  however,  is  that  the  case  burns  uniformly  along  its  length.  Erosive  effects  could 
increase  the  rate  of  surface  regression  at  the  case  mouth  (i.e.,  projectile  end  of  case), 
however,  tests  with  inert  cases  show  no  such  preferential  erosion..  On  the  other  hand,  the 
pressure,  which  in  the  code  is  independent  of  distance,  in  reality  is  higher  at  the  breech  and 
this  fact  tends  to  increase  the  rate  of  burning  at  the  breech.  To  some  unknown  extent,  these 
effects  compensate  for  each  other. 

6t  Combustion  of  the  fragments  is  not  considered.  Presumably,  combustion  would  increase 
the  rate  of  consumption  and  heating  of  the  particle,  so  that  here  an  upper  bound  to  the 
survival  time  of  the  fragment  is  computed. 


3.  Conservation  Equation  Governing  the  Fragment  Temperature 

If  r  is  the  time  after  injection  of  the  fragment  into  the  flow,  (other  symbols  are  identified  in 
the  List  of  Symbols) 

pVc^f  =  AhiTj-Tl+eoACl^-T4)  U) 


Neglecting  radiation  heat  transfer  (for  the  moment),  and  assuming  T ^  is  constant, 


pVc^f  =  Ah{Ttf-T) 
dx 


(2) 


Changing  variables  to  ©  =  -  T  and  considering  h  constant, 

—  =  —  q  -  —  Q 

dx  p  Vc  p  dc 

Assuming  T  =  Tq  at  r  =  0  ( T0  =  300  K  throughout  this  report.) 

— T 

Tgf-T  =  {T^-Tje 


(3) 


(4) 


2 


7W  -  V(Vro)«  ** 


(5) 


This  equation  will  govern  the  temperature  rise  of  the  fragment  up  to  the  melting  point, 
requiring  a  time  Atj  given  by 


Ax, 


=  Pg£|n (  Ij£12l 
2A 


(6) 


at  the  end  of  which  time  the  temperature  will  remain  constant  until  the  heat  of  fusion  is 
delivered  to  the  fragment.  This  requires  a  period  of  time  Arm  given  by 


At 


m 


2 h(T„-TJ 


(V 


When  melting  is  complete,  the  melted  mass  (assuming  it  remains  intact,  has  the  same  shape, 
density,  and  heat  capacity)  resumes  its  temperature  rise  according  to 


2 h  t 

no  ■  v<Vr> 


(8) 


for  a  time  interval  Ar2  given  by 


At2  =  ^Inl 
2  2 h 


\TgTTb ) 


(9) 


When  the  boiling  point  is  reached,  the  temperature  of  the  molten  fragment  again  remains 
constant  for  a  time  Arvap  until  the  heat  of  vaporization  is  delivered. 


At 


yap 


prfA//^ 

2  MVr>) 


(10) 


At  the  end  of  this  period  the  fragment  is  completely  vaporized. 


4.  Rate  of  Heat  Transfer  to  the  Fragment 

The  heart  of  the  fragment  heating  problem  is  in  the  computation  of  the  heat  transfer 
coefficient.  This  necessitates  an  empirical  correlation  function,  given  in  the  next  section,  and 
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a  calculation  of  the  relative  velocity  between  the  accelerating  fragment  ana  the  gas  flow, 
which  is  itself  a  function  of  time  and  position. 


4.1  Heat  Transfer  Coefficients 

Case  1:  Static  Gas 

Melvin1  determined  a  heat  transfer  coefficient  for  a  metal  sphere  in  a  static  gas.  For  a 
sphere  of  radius  a 


h 


(11) 


Defining  an  equivalent  radius  of  the  fragment  wafer  as  the  radius  of  a  sphere  of  same 
surface  area  as  the  fragment, 


a 


N  4u 


(12) 


Case  2:  Moving  Gas 

The  heat  transfer  coefficient  used  in  the  NOVA  interior-ballistic  code  is  due  to  Gelperin  and 
Ainstein2.  This  correlation  is  for  spheres.  Its  use  for  non-spherical  particles  of  equivalent 
area,  though  common,  has  not  been  tested. 

Nu  =  1+0.2 Re0t7Pr033  (13> 


where 


Nu 


(14) 


Re  • 


2P,Vfl 


(15) 


where  v  is  the  velocity  of  the  flow  relative  to  the  fragment. 


Pr  s 


Vi 

X 


(16) 


Eqs.  13  and  14  reduce  to  Eq.  11  in  the  limit  of  zero  flow  velocity. 


4.2  Computation  of  Flow  Conditions  at  the  Fragment 


The  fragment  location  in  the  cartridge  case  most  likely  to  result  in  the  fragment’:  survival 
is  the  outside  surface  of  the  case  in  contact  with  the  chamber  wall.  All  calculations  will  be 
based  on  this  radial  location  though  different  axial  locations  will  be  considered.  Thus  the 
fragment  is  considered  to  be  "injected"  into  the  interior-ballistic  flow  at  the  time  of  burnout 
of  the  combustible  case.  In  order  to  compute  the  temperature  history  of  the  fragment,  one 
must  estimate  the  velocity  of  the  combustion  gases  relative  to  that  of  the  fragment,  (v^--  vy). 
In  order  to  determine  this  velocity,  one  must  compute  the  trajectory  of  the  fragment  along 
the  flow  streamline.  This  calculation  proceeds  as  follows.  (Fig.  1  may  be  helpful  in 
visualizing  the  variables  at  the  time  of  injection,  /, ,  and  at  the  time  of  projectile  exit,  tc .) 


4.3  Equation  of  Motion  of  the  Fragment 


The  velocity  of  the  fragment,  Vy ,  with  mass  m  is  obtained  by  solving  its  equation  of  motion: 


(17) 


where  is  the  drag  coefficient  of  the  fragment ,  -v/)2/2  is  the  free-stream  dynamic 

pressure  on  the  fragment,  and  As  is  the  sectional  fragment  area  presented  to  the  flow.  For 
constant  v ^  the  fragment  equation  of  motion  can  be  solved  analytically  to  give 


v/t) 


1+Pvr 


(18) 


where  /?is  defined  by 


p  =  CpP  gAs 
2m 


(19) 


Note  that  Cq  has  been  assumed  constant  in  Eq.  17  but  that  it  often  is  a  function  of  the 
Reynolds  number,  particularly  at  the  transition  from  laminar  to  turbulent  flow.  The 
maximum  drag  on  the  fragment  occurs  when  the  normal  to  the  fragment  surface  is  oriented 
in  the  direction  of  the  flow.  The  value  of  the  drag  coefficient  under  this  condition  is 
approximately  1.  In  orientations  where  the  surface  normal  is  perpendicular  to  the  flow,  the 
drag  will  assume  a  very  low  value  due  to  the  small  edge  area  presented.  One  could  average 
the  drag  over  all  possible  orientations  under  the  assumption  that  the  particle  is  tumbling; 
however,  for  purposes  of  this  limited  analysis,  it  is  deemed  preferable  to  compute  bounds 
on  the  behavior  to  avoid  making  unwarranted  assumptions.  The  maximum  value  of  /?  is 
therefore  computed  using  a  Cp  =  1  and  As  =  s 2,  whereas  the  minimum  value  of  fi assumes 
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Cjy  =  1  and  As  =  sd.  Although  the  minimum  value  of  Cp  would  be  less  than  1,  the  very 
small  edge  area  should  provide  an  effective  lower  bound  on  0. 


Integrating  from  Eq.  18  over  r,  one  gets  for  the  position  of  the  fragment  Xj 


*/(*) 


=  w 


ln(l+pyx) 

P 


(20) 


The  approximation  that  y  is  constant  in  the  fragment  equation  of  motion  will  now  be 
examined.  The  velocity  of  the  combustion  gases  (more  precisely,  the  two-phase  combusting 
fluid),  Vg  ,  at  any  time  and  position  behind  the  projectile  is  given  by 


vg(x,t)  = 


xvp(t) 


(21) 


where  x  is  the  distance  from  the  breech  block,  vp(t)  is  the  projectile  velocity  at  time  t,  Lc  is 
the  length  of  the  chamber,  and  £p(t)  is  the  distance  traveled  by  the  projectile  at  time  t.  This 
linear  velocity  profile  is  consistent  with  the  assumptions  of  the  computer  code3  (IBHVG2) 
used  to  characterize  the  interior  ballistic  conditions.  The  time  r  is  referenced  to  the  time 
of  combustible-case  burnout,  i.e.,  the  fragment  is  injected  at  r  =  0,  or  f,-.  The  motion  of  a 
fluid  particle  originating  at  the  position  and  time  of  injection,  is  obtained  by  integrating 
Eq.  21  as  follows: 


Jdx  m  }  vfr'+t)  dx,  (22) 

*,  x  o  VMT'+r«) 


The  integrand  on  the  right  side  of  this  equation  is  shown  in  Fig.  2  plotted  as  a  function  of 
gun  time  t.  Case  burnout  occurs  at  f,  =  3.85  ms,  so  the  integrand  is  well  approximated  as 
a  linear  function  for  times  longer  than  3.85  ms.  Thus,  for  r  >  0,  based  on  a  simple  two- 
point  fit, 


— V-2—  =  587-1. 12x105t  (s'1)  (23) 

L  +1 

C  *• p 


Using  Eq.  23  in  Eq.  22,  the  position  of  the  fluid  particle  at  time  r  is  given  by 

JC,(t)  =  VM7t-5.6x10V 


(24) 
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Eqs.  18  and  20  were  derived  on  the  assumption  of  constant  v ^  We  now  discuss  how  to 
choose  a  suitable  value  for  purposes  of  this  analysis.  Since,  at  a  given  time,  the  gas  velocity 
is  a  linear  function  of  position,  a  lower  bound  on  the  velocity  of  the  gas  in  contact  with  the 
fragment  can  be  found  by  using  x,  in  Eq.  21.  Since  the  fragment  can  never  completely  catch 
up  with  the  fluid  particle  which  passes  it  at  the  instant  of  injection,  an  upper  limit  on  can 
be  constructed  by  using  Eq.  24  in  Eq.  21.  Thus,  upper  and  lower  bounds  on  v^as  a  function 
of  time  can  be  readily  computed  and  are  shown  in  Figs.  3  -  5  for  three  axial  injection 
positions.  The  heat  transfer  to  the  fragment  is  very  sensitive  to  injection  position  near  the 
breech  face  because  IBHVG2  assumes  that  the  velocity  is  zero  at  the  breech  wall.  In  reality 
even  minor  turbulence  will  sweep  the  fragment  from  the  breech  wall.  To  reflect  this 
behavior,  we  have  arbitrarily  taken  1  cm  as  the  minimum  distance  from  the  breech  to  be 
considered  here.  Calculations  assuming  that  is  constant  at  the  injection  value  v^j ,  i.e., 
v„  at  x,-  and  r  =  0,  can  be  used  to  calculate  (using  Eq.  20)  xy  at  projectile  exit  time,  rc  ,  and 
this  distance  can  in  turn  be  used  in  Eq.  21  to  estimate  v_^  at  projectile  exit  time.  We  label 
this  value,  v^c  ,  and  print  out  the  value  for  every  run.  In  the  cases  thus  far  examined  v^c 
is  within  30  %  of  in  almost  all  cases  and  always  within  37  %.  The  accuracy  of  tne 
computations  could,  of  course,  be  improved  by  numerically  integrating  the  coupled,  non¬ 
linear  differential  equations  describing  the  system,  but  in  our  judgement  this  error  is  by  far 
overshadowed  by  the  uncertainty  in  the  drag  on  the  fragment. 


4.4  Bounds  on  the  Heat  Transfer 

In  order  to  compute  the  heat  transfer  from  the  flow  to  the  fragment  after  injection,  the 
velocity  of  the  gas  in  contact  with  the  fragment  relative  to  the  velocity  of  the  fragment  is 
required,  i.e.,  v  in  Eq.  15,  where 

v  =  v  ,-v,  (25) 


We  will  compute  the  average  value  of  v  in  the  interval  (0,rc)  as  follows: 

1  T‘  1 

<V>  =  <Vgf-Vf>  =  <V>-<V,>  =  ±fVgfdX-±fvfdT 

0  Te0 


(26) 


<v> 


(27) 


vg  in  Eq.  26  is  taken  as  to  be  consistent  with  the  assumptions  of  the  fragment-trajectory 
calculation.  In  those  cases  where  increases  appreciably  over  the  time  interval  (0,r)  Vy 
should  also  increase,  tending  to  minimize  changes  in  v. 
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Max  and  Min  Gas  Velocities  after  Injection 
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Fig.  3.  Upper  and  lower  bounds  on  after  injection  near  breech 
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Max  and  Min  Gas  Velocities  after  Injection 
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Fig.  5.  Upper  and  lower  bounds  on  after  injection  at  case  mouth 


Since  the  derivations  of  Eqs.  5  -  10  assumed  constant  Tg  ,  we  may  maximize  the  heat 
transfer  by  choosing  as  Tg  the  maximum  value  of  the  flow  temperature  Tg  in  the  interval  (0, 
rc),  using  the  minimum  drag  constant  0 for  Xj  in  Eq.  20.  Similarly,  the  minimum  value  of  Tg 
may  be  used  for  Tg ,  coupled  with  the  maximum  0,  to  obtain  the  minimum  heat  transfer. 
The  arithmetic  average  of  these  two  extreme  temperatures  is  used  in  conjunction  with  the 
arithmetic  average  value  of  0  in  the  "average"  temperature  histories  shown  in  the  figures. 
It  should  be  pointed  out  that  use  of  the  average  relative  velocity  in  Eq.  26,  in  a  strict  sense, 
invalidates  the  temperature  history  details  computed  by  Eq.  2.  In  reality,  the  relative 
velocity  at  early  times  is  much  higher  than  <v>  and  at  late  times  much  lower.  However, 
the  computed  value  of  fragment  temperature  at  projectile  exit  time  should  be  valid,  since 
the  <vy>  in  Eq.  27  is  the  average  value  which  properly  takes  into  account  the  non-linear 
time  dependence  of  vy  given  by  Eq.  18. 


5.  Radiative  Heat  Transfer  vs.  Convective  Heat  Transfer 

Radiative  heat  transfer  can  be  shown  to  be  negligible  compared  to  that  by  convection  by 
examining  the  case  of  minimum  convective  transfer,  i.e.,  at  an  injection  position  near  the 
breech.  In  Fig.  6  the  fragment  thermal  history  determined  by  Eq.  2  is  used  to  compute  the 
relative  magnitude  of  these  two  heat  transfer  mechanisms.  The  radiative  transfer  is  dwarfed 
by  convection  even  under  the  minimum  convection  conditions.  This  justifies  Eq.  2  as  an 
excellent  approximation. 


6.  Results  and  Discussion 

All  of  the  results  thus  far  computed  pertain  to  the  M829  round.  Figs.  7  -  12  show  the 
temperature  histories  of  a  square  fragment  0.3  cm  on  a  side  by  0.03  cm  thick.  Figs.  7  and 
8  assume  the  fragment  (aluminum  and  steel,  respectively)  is  located  initially  at  the  case 
mouth.  Figs.  9  and  10  assume  the  fragment  is  initially  at  the  midlength  of  the  case,  and  Figs. 
11  and  12  take  the  initial  fragment  location  as  1  cm  from  the  breech  wall.  Clearly,  the 
assumptions  relating  to  the  drag  on  the  fragment  have  a  very  strong  effect  on  the  computed 
results,  particularly  at  the  case  mouth.  If  the  drag  is  high,  then  the  fragment  quickly 
accommodates  to  the  velocity  of  the  flow  and  the  convective  heat  transfer  is  minimized.  The 
volumetric  heat  capacity  of  steel  is  about  28  %  higher  than  aluminum  over  the  temperature 
range  considered.  Therefore,  for  equal  heating  rates  the  temperature  of  the  steel  is  lower 
than  that  of  the  aluminum.  However,  the  mass  of  the  steel  fragment  is  higher,  leading  to 
a  lower  acceleration,  and  therefore  higher  relative  velocity  between  the  fragment  and  the 
flow.  This  tends  to  increase  the  heat  transfer  to  the  sreel  fragment,  in  opposition  to  the 
heat-capacity  effect.  For  the  most  part,  the  increased  heat-transfer  effect  dominates,  leading 
to  higher  temperatures  for  the  steel  fragment.  Exceptions  to  this  pattern  are  often  found 
near  the  breech,  however  (e.g.,  Figs.  11  and  12).  Here  the  fragment  velocity  is  a  small 
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Fig.  6.  Comparison  of  radiative  to  conductive  heat  transfer 
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Fig.  7.  Temperature  history  of  aluminum  fragment  (s  =  .3  cm,  d  =  .03  cm)  originating  at 
the  case  mouth 
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Fig.  8.  Temperature  history  of  steel  fragment  (s  =  .3  cm,  d  =  .03  cm)  originating  at  the 
case  mouth 
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Fig.  9.  Temperature  history  of  aluminum  fragment  (s  =  .3  cm,  d  =  .03  cm)  originating  at 
case  midlength 


M829  :  Fragment  at  Case  Midlength 
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Fig.  10.  Temperature  history  of  steel  fragment  (s  =  .3  cm,  d  =  .03  cm)  originating  at  case 
midlength 


M829  :  Fragment  near  Breech 


[  x  ]  ajn^Bjadojai  }uaw6ejj 


19 


Fig.  11.  Temperature  history  of  aluminum  fragment  (s  =  .3  cm,  d  =  .03  cm)  originating 
near  the  breech 


M829  :  Fragment  near  Breech 
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Fig.  12.  Temperature  history  of  steel  fragment  (s  =  .3  cm,  d  =  .03  cm)  originating  near  the 
breech 


fraction  of  the  flow  velocity  for  both  aluminum  and  steel,  minimizing  the  heat  transfer  effect 
and  emphasizing  the  heat-capacity  effect. 

The  effect  of  both  increasing  and  decreasing  the  fragment  side  dimension  by  a  factor  of  2 
is  shown  in  Figs.  13  -  20.  In  general,  the  highest  fragment  temperatures  are  associated  with 
case-mouth  origins,  but  these  fragments  also  achieve  the  highest  forward  velocities  and 
therefore  are  likely  to  leave  the  gun  tube  after  shot  ejection.  It  is  interesting  to  note, 
however,  that  in  this  series  of  runs  the  fragment  position  at  the  time  of  projectile  exit  never 
reaches  more  than  25  %  of  the  combined  barrel/chamber  length  (531  cm)  because  of  the 
late  time  of  injection  during  the  IB  cycle. 

Decreasing  the  thickness  of  the  0.3  cm  side  fragment  to  0.02  cm  changes  the  thermal 
histories  as  shown  in  Figs.  21  -  26. 

A  set  of  fragment  dimensions  resulting  in  complete  vaporization  is  found  in  Fig.  27.  The 
fragment  thickness  here  is  0.0025  cm  (0.001  in.)  and  its  side  dimension  is  0.05  cm  (0.02  in.). 
This  size  is  very  close  to  the  threshold  of  complete  vaporization.  (The  dotted  line  ends  at 
the  time  of  complete  conversion  to  vapor.)  It  clearly  takes  a  relatively  long  time  to  supply 
the  large  (2720  cal/g)  heat  of  vaporization.  Under  "average"  assumptions  in  this  same  figure, 
the  fragment  has  not  quite  begun  to  vaporize.  Under  the  same  conditions,  a  steel  fragment 
(Fig.  28)  stabilizes  at  the  maximum  flow  temperatures  in  the  molten  state.  Figs.  29  and  30 
show  the  same  size  fragment  originating  near  the  breech.  Note  that  velocities  of  these 
fragments,  both  aluminum  and  steel,  are  low  but  their  temperatures  are  quite  high  and 
therefore  may  pose  a  serious  ignition  threat. 


7.  Conclusions 

The  problem  of  a  metallic  fragment,  trapped  in  a  combustible  cartridge  case,  and  heated  by 
the  combustion  gases  during  the  interior-ballistic  cycle  has  been  posed  mathematically  in  this 
report.  In  order  to  compute  the  convective  heat  transfer  from  the  combustion  gases,  one 
must  compute  the  trajectory  of  the  fragment  along  the  flow  streamline,  then  couple  this 
information  to  tiie  equation  governing  the  internal  energy  of  the  fragment  An  approximate 
solution  to  this  problem  has  been  constructed  using  analytic  solutions  valid  under  restricted 
conditions.  A  numerical  solution  to  the  set  of  coupled  non-linear  differential  equations 
describing  this  system  would  eliminate  the  inaccuracies  of  the  analytic  approach  used; 
however,  in  our  judgement,  the  largest  uncertainties  in  these  calculations  arise  from 
uncertainties  in  the  description  of  the  drag  of  a  wafer-like  particle  tumbling  in  the  flow.  In 
the  calculations  presented  here,  this  uncertainty  is  treated  by  computing  bounds  on  the 
thermal  history  of  the  fragment  arising  from  bounds  on  the  fragment  drag.  It  could  well  be 
that  the  chaotic  motion  of  fragments  in  the  interior-ballistic  flow  would  be  characterized  by 
large  fluctuations  about  the  mean  behavior.  If  this  is  the  case,  placing  bounds  on  the 
behavior  is  the  most  meaningful  approach  from  an  engineering-design  perspective. 
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M029  :  Fragment  at  Case  Mouth 


cd 

o 

CO 


CD 

CD 

OJ 


CD 

CD 


CD 

CD 

CD 


CD 

CD 

CD 

m 


CD  CD  CD 

CD  CD  CD 

CD  CD  CD 

OJ  *c-l 


c n 


[  )\  ]  ajn^Bjadiuai  }uaui6ejj 


22 


Fig.  13.  Temperature  history  of  aluminum  fragment  (s  =  .6  cm,  d  =  .03  cm)  originating  at 
case  mouth 


M829  :  Fragment  at  Case  Mouth 
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Fig.  14.  Temperature  h'story  of  steel  fragment  (s  =  .6  mm,  d  =  .03  mm)  originating  at  the 
case  mouth 
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Fig.  15.  Temperature  history  of  aluminum  fragment  (s  =  .6  mm,  d  =  .03  mm)  originating 
near  the  breech 


M829  :  Fragment  near  Breech 
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Fig.  16.  Temperature  history  of  steel  fragment  (s  =  .6  cm,  d  =  .03  cm)  originating  near  the 
breech 


M829  :  Fragment  at  Case  Mouth 
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Fig.  17.  Temperature  history  of  aluminum  fragment  (s  =  .15  cm,  d  =  .03  cm)  originating 
at  the  case  mouth 


M829  :  Fragment  at  Case  Mouth 
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Time  after  Injection  [  ms  ]  0B/ 

Fig.  18.  Temperature  history  of  steel  fragment  (s  =  .15  cm,  d  =  .03  cm)  originating  at  the 
case  mouth 
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Fig.  19.  Temperature  history  of  aluminum  fragment  (s  =  .15  cm,  d  =  .03  cm)  originating 
near  the  breech 


M829  :  Fragment  near  Breech 
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Fig.  20.  Temperature  history  of  steel  fragment  (s  =  .15  cm,  d  =  .03  cm)  originating  near 
the  breech 


[  x  ]  9jn}ej3dujai  ;u9qj6bjj 


30 


Fig.  21.  Temperature  history  of  aluminum  fragment  (s  =  .3  cm,  d  =  .02  cm)  originating  at 
case  mouth 
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Fig.  22.  Temperature  history  of  steel  fragment  (s  =  .3  mm,  d  =  .02  mm)  originating  at  the 
case  mouth 


M829  :  Fragment  at  Case  Midlength 
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Fig.  23.  Temperature  history  of  aluminum  fragment  (s  =  .3  mm,  d  =  .02  mm)  originating 
at  case  midlength 


M829  :  Fragment  at  Case  Midlength 
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ig.  24.  Temperature  history  of  steel  fragment  (s  =  .3  cm,  d  =  .02  cm)  originating  at  case 
lidlength 


MB29  :  Fragment  near  Breech 
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Fig.  25.  Temperature  history  of  aluminum  fragment  (s  =  .3  cm,  d  =  .02  cm)  originating 
near  the  breech 


M829  :  Fragment  near  Breech 
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Fig.  26.  Temperature  history  of  steel  fragment  (s  =  .3  cm,  d  =  .02  cm)  originating  near  the 
breech 
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Fig.  27.  Temperature  history  of  aluminum  fragment  (s  =  .05  cm,  d  =  .0025  cm)  originating 
at  the  case  mouth 
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Fig.  28.  Temperature  history  of  steel  fragment  (s  =  .05  cm,  d  =  .0025  cm)  originating  at 
the  case  mouth 
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Fig.  29.  Temperature  history  of  aluminum  fragment  (s  =  .05  cm,  d  =  .0025  cm)  originating 
at  the  breech 


M829  :  Fragment  near  Breech 
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Fig.  30.  Temperature  history  of  steel  fragment  (s  =  .05  cm,  d  -  .0025  cm)  originating  at 
the  breech 


The  calculations  reported  here  are  intended  to  examine  the  conditions  under  which  the 
fragments  are  most  likely  to  survive  vaporization  and,  hence,  the  fragments  are  assumed  to 
be  imbedded  at  the  outer  radius  of  the  case,  though  different  axial  locations  are  considered. 
Results  of  these  calculations  show  that,  for  this  outer-radius  position,  only  very  small 
fragments,  near  the  limits  of  unaided  visual  inspection,  will  undergo  complete  vaporization. 
Particles  imbedded  near  the  case  mouth  will  achieve  the  highest  temperatures  due  to  the 
high  gas  velocity  there;  however,  these  particles  will  also  have  ihc  highest  forward  velocity 
at  the  time  of  projectile  exit  and  therefore  are  most  likely  to  be  ejected  from  the  barrel. 
Fragments  imbedded  near  the  breech  are  the  ones  most  likely  to  be  retained  in  the  chamber 
after  functioning  of  the  propelling  charge  and  these  fragments  will  present  a  spectrum  of 
ignition  threats  from  smaller  particles  with  relatively  low  internal  energy  at  very  high 
temperatures  (>2000  K)  to  larger  particles  with  relatively  high  internal  energy  but  moderate 
temperatures  ( ~  500  K).  Ignition  probabilities  for  particles  up  to  1275  K  in  contact  with 
combustible-case  materials  could  be  measured  using  the  Hot  Fragment  Conductive  Ignition 
Test4,5  (HFCIT)  developed  at  BRL  under  the  LOVA  program.  The  test  shows  an  inverse 
relationship  between  fragment  mass  and  its  temperature  required  to  ignite  a  given  energetic 
material.  The  smallest  fragment  used  in  the  HFCIT  is  a  130  mg  steel  ball  bearing. 
Unpublished  records  at  BRL  from  approximately  the  1981  -  1983  period  show  that  this 
"fragment"  ignited  U.S.  combustible-case  material  at  a  temperature  of  850  K.  A  wafer  shape 
likely  would  require  a  higher  temperature  due  to  its  higher  surface-to-volume  ratio.  New 
techniques  would  have  to  be  developed  to  measure  the  ignition  potential  of  fragments  hotter 
than  1275  K. 
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APPENDIX  1:  Interior  Ballistic  Calculations  for  the  M829 
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APPENDIX  1:  Interior  Ballistic  Calculations  for  the  M829 


Results  of  IBHVG2  calculations3  for  the  M829  round,  performed  at  BRL,  are  shown  in  Figs. 
A1  -  A6.  This  is  a  one-dimensional  code  that  assumes  spatially  uniform  but  time-varying 
pressures  and  temperatures  in  the  gun  chamber. 
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Fig.  A-2.  M829 
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Fig.  A-6.  M829 
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APPENDIX  2:  Numerical  Parameter  Values  Used  in  Computations 
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APPENDIX  2;  Numerical  Parameter  Values  Used  in  Computations 

Interior  Ballistic  Conditions3: 

P  =  3670  atm  mean  pressure  at  combustible-case  burnout 
T  =  2610  K  mean  temperature  at  combustible-case  burnout 
Vp  =  3525  ft/s  =  1.07E5  cm/s  projectile  velocity  at  case  burnout 
r  =  1100  atm  mean  pressure  at  projectile  exit 
T  =  2090  K  mean  temperature  at  projectile  exit 
Vp  =  5481  ft/s  =  1.67E5  cm/s  projectile  velocity  at  projectile  exit 
tc  =  2.43  ms  from  time  of  case  burnout  to  projectile  exit 

Gas  Properties  at  Particular  Conditions6: 

X  =  .0056  cal/cm-s-K  for  JA2  @  2610  K  &  3670  atm 
X  =  .0046  cal/cm-s-K  for  JA2  @  2090  K  &  1100  atm 
H  =  .0080  g/s-cm  for  JA2  @  2610  K  &  3670  atm 
/i  =  .0068  g/s-cm  for  JA2  @  2090  K  &  1 100  atm 
Pg  =  .427  g/cc  for  JA2  @  2610  K  &  3670  atm 
Pg  =  .160  g/cc  for  JA2  @  2090  K  &  1100  atm 

Average  Gas  Properties: 

X  =  .0051  cal/cm-s-K 
H  =  .0074  g/s-cm 
Pg  =  0.294  g/cc 

Fragment  Properties: 


Aluminum: 

p  =  2.71  g/cc  Alloy  1100  [Ref.  7] 

c(T)  =  0.185  +  1.18E-4  T  cal/g-K  in  range  300  -  950  K  [Ref.  8] 
c  =  0.340  cal/g-K  avg.  value  between  Tq  and  Tb,  extrapolating  previous  equation 
Tm  =  1220  F  =  660  C  =  933  K  melt.  temp.  [Ref.  7] 

Tb  =  3740  F  =  2060  C  =  2333  K  boiling  pt.  [Ref.  7] 

AH  fusion  =  90  cal/g  heat  of  fusion  [Ref.  8] 

AHvap  =  2720  cal/g  heat  of  vaporization  [Ref.  8] 
f  =  .05  @  ~  1000  K  for  cleaned  surface  [Ref.  8] 

Steel: 

p  -  7.85  g/cc  steel,  AISI  C1020  (hot  worked)  [Ref.  7] 

c  =  0.15  cal/g-K  avg.  value  between  T0  and  Tc  for  plain  carbon  steel  0.2  -  0.6  %  C  [Ref.  8] 
Tm  =  2750  F  =  1510  C  =  1783  K  melt.  temp,  steel  [Ref.  7] 

Tm  =  1400  -  1765  K  for  plain  carbon  steel  2.0  -  0.1  %  C  [Ref.  8] 

Affusion  =  65  cal/g  heat  of  fusion  for  iron  [Ref.  8] 

Constants: 

a  =  1.356E-12  cal/cm2-s-K4  Stefan-Boltzmann  const. 
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radius  of  sphere  of  surface  area  A 

fragment  heat  capacity  (cal/g-K) 

fragment  thickness  (cm) 

heat  transfer  coefficient  (cal/cm2-K-s) 

fragment  side  dimension  (assuming  square  shape) 

time,  referenced  to  beginning  of  interior-ballistic  cycle 

time  of  projectile  exit  from  barrel 

time  of  fragment  injection  into  interior-ballistic  flow 

fragment  velocity 

velocity  of  combustion  gases 

velocity  of  combustion  gases  at  the  fragment  position 
velocity  of  projectile 

position  in  chamber/tube  relative  to  breech  block 
location  of  fragment 

location  of  gas  particle  (combustion  fluid)  which  passes  xi  at 

location  of  fragment  at  injection 

fragment  surface  area 

fragment  sectional  area 

fragment  heat  of  fusion  (cal/g) 

fragment  heat  of  vaporization  (cal/g) 

length  of  chamber 

Nusselt  number 

Prandtl  number 

Reynolds  number 

temperature  of  fragment  (K) 

fragment  boiling  point  (K) 

temperature  of  the  gas  at  some  specified  position  in  the  chamber/bore 

temperature  of  the  gas  in  contact  with  the  fragment 

fragment  melting  point  (K) 

initial  temperature  of  fragment 

=  Ad/2,  fragment  volume 

fragment  emissivity 


=  TgTT 

thermal  conductivity  of  the  combustion  gases 
viscosity  of  combustion  gases 
distance  traveled  by  projectile 
fragment  density  (g/cc) 
density  of  combustion  gases 
Stefan-Boltzmann  constant 

time,  reference  to  instant  of  fragment  injection  into  flow 
time  of  projectile  exit 
average  of  enclosed  variable 
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